Role of the d — f Coulomb interaction in intermediate valence and Kondo systems: a 

numerical renormalization group study 
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Using numerical renormalization group method the temperature dependences of the magnetic 
susceptibility xC^) s-nd specific heat C{T) are obtained in the single-impurity Anderson model 
with inclusion of the d — f Coulomb interaction. It is shown that the exciton effects owing to 
this interaction can change considerably the dependence C(T) in comparison with the standard 
Anderson model at not too low temperatures, whereas the dependence x(T) remains universal. The 
renormalization of the effective hybridization parameter and /-level position, which is connected with 
the d — / interaction, is calculated, a satisfactory agreement with the Hartree-Fock approximation 
being derived. 
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There is an interesting class of the 4/-electron com- 
pounds demonstrating intermediate valence (IV) of rare 
earth elements (usually between 2-1- and 3-I-) in a num- 
ber of properties, e.g., in the lattice constants (which are 
intermediate between those for isostructural compounds 
with di- and trivalent ions), core- level spectra (which 
are mixtures of the spectra for di- and trivalent ions 
with comparable weights) , and many others Q, S S ■ 
Heavy fermion (HF) compounds "F] form another impor- 
tant class of the /-electron systems with anomalous prop- 
erties. For the HF metals it is commonly accepted now 
that they are the Kondo lattices, which means that the 
small energy scale in the electron properties is the Kondo 
temperature T/f, i.e. the width of the Kondo resonance 
owing to spin-dependent scattering of conduction elec- 
trons by /-electron centers 0. As for the IV compounds, 
they are frequently considered also as the Kondo lattices, 
but just with higher Tk (see, e.g., Ref|3). 

Actually, such a consideration is not quite accurate 
since, besides the spin ("Kondo") fluctuations, valence 
or charge fluctuations should be also treated in such sys- 
tems. They are determined in part by the Coulomb re- 
pulson G between conduction and localized electrons (the 
Falicov-Kimball interaction 7] ) . Taking into account the 
d — f interaction together with the hybridization pro- 
cesses it is possible to describe the IV state as a kind of 
exciton condensation Recently, the method of first- 

principle calculations of the parameter G has been pro- 
posed, and it was demonstrated that an account of this 
interaction is necessary to describe properly the equation 
of state for IV phase of Yb under pressure ^3 ■ 

At present, the usual Kondo effect is theoretically stud- 
ied thoroughly within the s — d exchange (Kondo) and 
Anderson models. Moreover, in the one-impurity situa- 
tion the exact numerical (renormalization group) [ill 
and analytical Bethe-ansatz jl^ J4J solution of this prob- 
lem is obtained. Universal curves describing the behav- 
ior of thermodynamic properties were obtained for the 



Kondo 0, 0| and intermediate valence ^3 regimes, 
which permit a detailed comparison with experimental 
data on anomalous /-systems. At the same time, in the 
presence of both the s — d exchange and Coulomb inter- 
action such a detailed information is absent. 

Formally, the charge fluctuations can be also described 
in terms of a pseudo-Kondo effect, the states with (with- 
out) /-hole being considered as pseudospin-up (down) 
states, respectively It is the degeneracy of quan- 

tum states for a scattering center which is important for 
the formation of the Kondo resonance the IV 

case the divalent and trivalent states are degenerate by 
definition, thus this analogy is not surprising. There- 
fore it is natural to consider the Kondo phenomenon for 
the IV compounds taking into account both spin and 
charge fluctuations, or, equivalently, both the "Kondo" 
and exciton ("Falicov-Kimball") effects. Since there is no 
clear demarcation between the IV and Kondo systems, it 
can be supposed that the exciton effects are relevant also 
for the latter case. Recent analysis of the interplay of 
the true Kondo and pseudo-Kondo (exciton) effects 0| 
by the "poor-man scaling" approach 0, [ifl] has demon- 
strated an essential modification of the low-energy (in- 
frared) behavior in comparison with pure cases of the 
Anderson model and Falicov-Kimball ("resonant level") 
models. However, this approach can give only a quali- 
tative insight in the properties of the system. Here we 
investigate the effects of this interplay by appl ying nu- 
merical renormalization group (NRG) approach |llL Il2l | . 

We proceed with the Hamiltonian of the asymmetric 
Anderson model with inclusion of the Falicov-Kimball 
interaction (on-site d — f Coulomb repulsion G) , 
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where the on-site f — f Coulomb interaction U is put to 
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infinity, so that the doubly occupied states are forbidden; 
fL — the Hubbard operator {\ia) and \iO) 

are single-occupied and empty site states); we neglect 
for simplicity k-dependence of the hybridization matrix 
element V. 

We use the standard NRG method for the Anderson 
model |l2l| with some important modifications. At each 
NRG step one obtains a finite-resolution spectrum which 
is truncated owing to neglect of high-energy states 
Thus we have a sequence of truncated energy spectra, 
the level resolution decreasing with increasing iteration 
step. An automatic choose of an optimal temperature 
for each NRG step is a main distinctive feature of our 
calculations. Indeed, we cannot perform calculations of 
thermodynamic averages at too low temperatures since 
the discreteness of energy level leads to an uncontrollable 
error. On the other hand, at sufficiently high tempera- 
tures the high-energy states neglected can give an ap- 
preciable contribution to the partition function, which is 
proportional to the factor of exp (— S/fcsT). Therefore 
we estimate the contribution of the upper states and in- 
crease the temperature to make their contribution to be 
equal to the error chosen. 

We have used a rectangular conduction electron den- 
sity of states with the half-bandwidth of Z) — 1. The 
magnetic susceptibility, specific heat, and impurity level 
occupation number Uf (valence) were calculated. Com- 
putational results are shown in the figures. Figs. ^ 
and 121 demonstrate a crossover from a two-maximum to 
one-maximum temperature behavior of specific heat. It 
should be noted that such a crossover takes place also in 
the standard Anderson model with changing Ef plj . 

TABLE I: The dependences of the Kondo temperature Tk 
and impurity occupation number Uf on G for V = 0.1; the 
Hartree-Fock values V^^ and E^^ are defined by Eq. |(2J. 
The quantities £;'°"°' and v'-°=°^ are discussed in the text. 
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FIG. 1: Effective Curie constant ksTxiT) / (gfiB)^ and spe- 
cific heat C{T)/kB for Ef = -0.14, V = 0.1. Solid line corre- 
sponds to finite G, and the dotted line to the case G=Q with 
parameters Z?^*^""' and l/C^"' (see their values in Table EJ- 
Below Tk we have the universal Wilson curve. 
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FIG. 2: The dependences of n/(«) and Tk{o) on G. {Ef = 
-0.14,1/ = 0.1). 
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One can see that the temperature dependence of the 
magnetic susceptibihty is always similar to that in the 
Anderson model without d — f Coulomb interaction. At 
the same, the specific heat behavior can be considerably 
different, especially for sufficiently large G. 

Similar to RefEl the Kondo temperature Tk was de- 
termined from the condition kBTx{T) / [g^sY ^ 0.0701 
at T = Tk- The dependence of the ground-state occupa- 
tion number of /-level n / and the Kondo temperature on 
G are illustrated by Fig. 13 A more detailed information 
is presented in the Tabled 

There is an important question whether the effects of 
the d ~ f Coulomb interaction can be described just by 
the renormalization of the parameters usual Anderson 
Hamiltonian (without the Falicov interaction) or they 
can result in qualitatively new effects. To investigate this 
problem we defined the effective hybridization parameter 
Y{g=o) ^-^^ .(.jjg effective position of the /-level, £'|*^~°\ 
as the parameters of the standard Anderson model (with 
G = 0) that gives the same values of n/ (at zero temper- 
ature) and the Kondo temperature, as our Hamiltonian 
(1). A comparison of our computational results with 
those for the model with G — and with the effective 
parameters introduced above shows (Fig. ElandEJ that 
for the susceptibility the effects of G in the temperature 
dependence can be completely eliminated by the param- 
eter renormalization. At the same time, for the specific 
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FIG. 4: Effective Curie constant ksTxiT)/ (q^b)^ and spe- 
cific heat C{T)/kB for Ef^-0.3, V^O.l, G=0.3 (solid) and 
E 
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-0.0475, V'('^=°) =0.109, G=0 (dotted). For these pa- 
rameter sets we have n/=0.77. Insert shows the dependence 
y^^(r) according to Eq.©. 
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FIG. 3: The same data as in Fig. [T]for Ej = -0.25, V = 0.1. 0-25- 
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FIG. 5: Effective Curie constant fcsTx(T)/(<7Ms)^ and spe- 
cific heat C{T)/kB for £;/=-0.5, ^=0.1, G=0.4 (solid) and 
-0.14, V'-'^=°'>=0.1, G=0 (dotted). For these param- 
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eter sets we have n/=0.96 at T—Q. Insert shows the depen- 
dence V"''{T). 



heat this is, generally speaking, possible for low enough 
temperatures, of order of Tk or below. This means that 
the Wilson ratio is not infiuenced by the d — f inter- 
action at r < Tk, but its "temperature dependence" 
at higher temperatures is different for the cases G = 
and G 7^ 0. Of course, it is not surprising that the d — f 
Coulomb (but not exchange) interaction is less important 
for the magnetic susceptibility (which is connected only 
with spin degrees of freedom) than for the specific heat 
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(which characterizes both spin and charge fluctuations). 

It is interesting to compare our renormahzed model pa- 
rameters with their values from the unrestricted Hartree- 
Fock approximation 0, 

Ef^ = Ef + G^(ct c.), V"^ = V- G{cif,). (2) 

A comparison of the parameters of the effective Ander- 
son model and of the Hartrcc-Fock values (2), which are 
presented in the Tabled shows that this approximation 
works well enough, at least for not too large d — / in- 
teraction (G < 0.25). Thus corresponding Coulomb cor- 
relation effects are not important. This justifies the im- 
plementation of the Hartree-Fock approximation into the 
first-principle electronic structure calculations in RefE 
The dependences V^^ (T) according to Eq. (j2l are shown 
in inserts in Figs.Q, One can see that in the Kondo 
regime a maximum occurs which is qualitatively similar 
to the result of the poor-man scaling consideration ^1§\ . 

To conclude, we have obtained an accurate NRG so- 
lution of the one-impurity Anderson model with inclu- 
sion of the Falicov-Kimball interaction (excitonic effects) . 
Some new features in comparison with the standard An- 
derson model (in particular, in the temperature depen- 
dence of specific heat) occur. A generalization of the re- 
sults to a lattice case would be of interest for the theory 
of the Kondo lattices and IV compounds. 
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